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We point out that the B — ^ pp data have seriously constrained the possibility of resolving the B — > 
7T7T puzzle from the large observed B° — > 7r°7r° branching ratio in the available theoretical approaches. 
The next-to-leading-order (NLO) contributions from the vertex corrections, the quark loops, and 
the magnetic penguin evaluated in the perturbative QCD (PQCD) approach have saturated the 
experimental upper bound of the B° — > p°p° branching ratio, and do not help. The NLO PQCD 
predictions for the B° — > p^ p and B ± — > p ± p° branching ratios are consistent with the data. The 
inclusion of the NLO jet function from the soft-collinear effective theory into the QCD-improved 
factorization approach, though enhancing the B° — > 7r°7r° branching ratio sufficiently, overshoots 
the bound of the B° — > p°p° branching ratio, and deteriorates the predictions for the B ± —* n°K ± 
<—j ' and B° — > ir^K ± direct CP asymmetries. 
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I. INTRODUCTION 



The observed direct CP asymmetries and branching ratios of the B — > ttK, tttt decays pj, 

o 
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= (4.9 ±0.4) x 10" 6 , 
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= (1.45 ±0.29) x 10~ 6 
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Oh! B{B° -> ttV) = (1.45 ±0.29) x 10~ 6 , (1) 

were regarded as puzzles, since they obviously contradict to the expected relations Acp{B a — > tt^K^) w 
A CP (B ± — > -k^K^ 1 ) and B(B° — > tt^tt^ » B(B° — > 7r°7r°). These puzzles have been analyzed in the per- 
turbative QCD (PQCD) approach 0,0 up to next-to-leading-order (NLO) accuracy recently Q, where the 
contributions from the vertex corrections, the quark loops, and the magnetic penguin were taken into account. 
It was found that the vertex corrections modify the color-suppressed tree contribution, such that the relative 
strong phase between the tree and penguin amplitudes involved in the B — > txK decays decreases. The predicted 
magnitude of the B^ — > direct CP asymmetry then becomes smaller, and matches the data in Eq. (JTJ. 

Though the B — > irK puzzle has been resolved, the B — > 7T7t puzzle remains, because the NLO color-suppressed 
tree amplitude does not increase the predicted B° — > 7r°7r° branching ratio sufficiently. 

A resolution to a puzzle usually demands an introduction of new mechanism. It is thus essential to investigate 
whether the proposed new mechanism deteriorates the consistency of theoretical results with other data. To 
make sure the above NLO effects are reasonable, we apply the same PQCD formalism to more two-body 
nonleptonic B meson decays, concentrating on the B —> pp branching ratios, which are also sensitive to the 
color-suppressed tree contribution. It will be shown that the NLO PQCD predictions are in agreement with 
the data of the B° — > p T p^ 1 and B^ — ► p^p branching ratios, and saturate the experimental upper bound of 
the B° -» p a p° branching ratio, B{B° -> p a p°) < 1.1 x lO" 6 0. Therefore, our resolution to the B — > ttK 
puzzle makes sense, and the B — > 7T7t puzzle is confirmed. The dramatic difference between the B — > 7T7t and pp 
data has been also noticed in [jj, which stimulates the proposal of a new isospin amplitude with I — 5/2. The 



possible new physics signals from the B — > 7T7t decays have been discussed in 
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It has been claimed that the B — * tttt puzzle is resolved in th e Q CD-improved factorization (QCDF) approach 
with an input from soft-collinear effective theory (SCET) 01 : the inclusion of the NLO jet function, one 
of the hard coefficients of SCETn, into the QCDF formula for the color-suppressed tree amplitude leads to 
enough enhancement of the B — > 7r°7r° branching ratio. Following the argument made above, we apply the 
same formalism ^(J to the B — > nK, pp decays as a check. It turns out that the effect of the NLO jet function 
deteriorates the QCDF results for the direct CP asymmetries in the B^ — > it^K^ and B a — > tt^K^ decays: 
the magnitude of the former increases, while that of the latter decreases, contrary to the tendency indicated 
by the data. This NLO effect also overshoots the upper bound of the B° — ► p°p° branching ratio very much. 
This observation is expected: the B° — > p°p° and B° — ► tt q tt° decays have the similar factorization formulas, so 
the branching ratio B(B° — > p°p°) ought to be larger than B(B° — > 7r°7r°) due to the meson decay constants 
f p >fir- Therefore, the B — > pp data have seriously constrained the possibility of resolving the B — > tttt puzzle 
in the available theoretical approaches. 

There exists an alternative phenomenological application of SCET where the jet function, charac- 

terized by the scale of 0(y / mbA), m& being the b quark mass and A a hadronic scale, is regarded as being 
incalculable. Its contribution, together with other nonpcrturbative parameters, such as the charming penguin, 
were then determined by the B — ► rnr data. That is, the color-suppressed tree amplitude can not be explained, 
but the data are used to fit for the phenomenological parameters in the theory. Predictions for the B — > irK, 
KK decays were then made based on the obtained parameters and partial SU(3) flavor symmetry |l2j |. Final- 
state interaction (FSI) is certainly a plausible resolution to the B — > irir puzzle, but the estimate of its effect 
is quite model-dependent. Even opposite conclusions were drawn sometimes. When including FSI either into 
naive factorization |l3| or into QCDF |l4j, the B° — ► 7r°7r° branching ratio was treated as an input in order 
to fix the involved free parameters. Hence, no resolution was really proposed. It has been found that FSI, 
evaluated in the Regge model, is insufficient to account for the observed B° — > 7r°7r° branching ratio [l^. We 
conclude that there is no satisfactory resolution in the literature: the available proposals are either data fitting, 
or can not survive the constraints from the B — ► ttK, pp data under the current theoretical development. 

In Sec. II we compute the branching ratios, the direct CP asymmetries, and the polarization fractions of the 
B — > pp decays using the NLO PQCD formalism. The branching ratios and the direct CP asymmetries of the 
B — > irK, 7T7T, pp decays are calculated in Sec. Ill by including the NLO jet function from SCETn into the 
QCDF formulas. Section IV is the discussion, where we comment on and compare the various analyses of the 
FSI effects in the B — > irK, irir decays. 



II. B -> pp IN NLO PQCD 



The NLO contributions from the vertex corrections, the quark loops, and the magnetic penguin to the 
B — > irK and hit decays have been calculated in the naive dimensional regularization (NDR) scheme in the 
PQCD approach and the results for the branching ratios and the direct CP asymmetries are quoted in 
Tables [I] and [HJ respectively. We have taken this chance to correct a minor numerical mistake in the vertex 
corrections for the B — > ttK decays, whose branching ratios become smaller by 2 ~ 4%. Note that a minus sign is 
missing for the q = t term in the expression for the quark-loop contributions in Eq. (27) of 0. Nevertheless, this 
typo has nothing to do with the numerical outcomes. Our observations are summarized below. The corrections 
from the quark loops and from the magnetic penguin come with opposite signs, and sum to about —10% of the 
leading-order (LO) penguin amplitudes. They mainly reduce the penguin-dominated B — > irK branching ratios, 
but have a minor influence on the tree-dominated B — > tttt branching ratios, and on the direct CP asymmetries. 
On the contrary, the vertex corrections do not change the branching ratios, except the B° — ► ir°ir° one. They 
modify only the direct CP asymmetries of the B ± — > ir°K ± , B° — ► tt°K°, and B° — ► 7r°7r° modes by increasing 
the color-suppressed tree amplitude C few times. The larger C', leading to the nearly vanishing direct CP 
asymmetry Acp(B ± — * tt q K ± ), resolves the B — * ttK puzzle within the standard model. 

The above observations can be easily understood as follows. The B° — > 7r =F i4T ± decays involve the color- 
allowed tree T' and the QCD penguin P' in the topological amplitude parametrization. The data of Acp(B° — > 
7r =F J ftT ± ) f=a —11.5% imply a sizable relative strong phase between T 1 and P' . The S ± — > tt^K^ 1 decays involve 
C and the electroweak penguin amplitude P' ew , in addition to T' and P' . If C is large enough, and more or less 
orthogonal to T', it may orient the sum T" + C roughly along with P' + P' ew - The smaller relative strong phase 
between T' + C and P' + P' ew then gives A C p{B ± -> ^K^) sa 0. We found in PQCD that the vertex corrections 
indeed modify C in this way. Because our analysis shows the sensitivity of C to the NLO corrections, it is 
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Mode Data [1] LP LO NLO wc +VC +QL +MP +NLQ 

30.1 34.2 24.1 23.6+" [l 

17.1 19.4 14.0 13.6+ 1 °^ ( ( + J 3] 

26.1 29.4 20.5 20.4 +16 1 (+11 ' 5) 



B ± ^ 


t^X 


24.1 ± 1.3 


17.0 


32.3 


B ± -» 




12.1 ±0.8 


10.2 


18.4 


S° — » 




18.9 ±0.7 


14.2 


27.7 


B° — > 


7^° 


11.5 ± 1.0 


5.7 


12.1 


S° — » 




4.9 ±0.4 


7.0 


6.8 


B ± ^. 




5.5 ±0.6 


3.5 


4.1 


B° — f 




1.45 ±0.29 


0.12 


0.27 



11.4 12.8 8.7 8>7 + e-o (+5.5) 



8.4 (- 6.7) 
6.0 (+ 5.5) 
3.4(- 3.1) 



6.6 6.9 6.7 6.5 + 6 ' 7(+27) 



4.0 4.1 4.1 4.0 +3 ' 4(+1 - 7) 
0.37 0.29 0.21 o.29 +0 ' 50(+0 13) 



3.8(- 1. 
3.4 (+ 1.' 
1.9 (- 1.2) 
).S0(+0J 
0.20 (-0.1 



TABLE I: Branching ratios from PQCD in the NDR scheme in units of 10 -6 . The label LOnlowc means the LO results 
with the NLO Wilson coefficients, and +VC, ±QL, ±MP, and ±NLO mean the inclusions of the vertex corrections, of the 
quark loops, of the magnetic penguin, and of all the above NLO corrections, respectively. The errors in the parentheses 
represent only the hadronic uncertainty |4j. 



Mode Data [1] LO LO NLO wc ±VC ±QL ±MP ±NLO 



B ± - n ± K° 


-2±4 


-1 


-1 


-1 





-1 


± (±0) 


B ± tt ^ 
B° -> - 
B° Tf K° 


4±4 
-10.8 ± 1.7 
2 ±13 


-8 
-12 

-2 


-6 
-8 



-2 
-9 
-7 


-5 
-6 



-8 
-10 



1+3 (+3) 
-6(-5) 
ln +7(+5) 
±U -8(-6) 
7+3 ( + 1) 
'-4-2) 


B° -> TT^ 


37 ± 10 


14 


19 


21 


16 


20 


1 „+20 (+ 7) 
10 -12 (- 6) 


B ± - T^TT 


1 ±6 

















0±0(±0) 


B° - ttV 


OO+40 
28 -39 


-4 


-34 


65 


-41 


-43 


fi „+35 (+ 9) 
UO -34 (-15) 



TABLE II: Direct CP asymmetries from PQCD in the NDR scheme in percentage. 



worthwhile to investigate the direct CP asymmetries of other charged B meson decays. The results will be 
published elsewhere. The color-suppressed tree amplitude C involved in the B —> irir decays, despite of being 
increased few times too by the vertex corrections, remains subleading with the ratio \C/T\ « 0.2, where T 
represents the color-allowed tree amplitude. This ratio is not enough to explain the observed B° — > 7r°7r° 
branching ratio as shown in Table |I] 4]. A much larger \C/T\ w 0.8 must be achieved in order to resolve the 
B — > 7T7T puzzle 0| . We mention that a different source for the large relative strong phase between C and T 
has been proposed in 0], which arises from charm- and top- mediated penguins. 



A. Helicity Amplitudes 



We examine whether the observations made in |4J are solid by applying the same NLO PQCD formalism to 
the B — > pp decays, which are also sensitive to the color-suppressed tree contribution. The B — > pp decays have 
been analyzed at LO in 0, 0] . The numerical results in the two references differ a bit due to the different 
choices of the characteristic hard scales, which can be considered as one of the sources of theoretical uncertainties 
(from higher-order corrections). The B — > Va(e2, ^2)^3 (£3, P3) decay rate is written as 

B a 

where P c — |Pa| = IP3I = ms/2 is the momentum of either of the vector mesons V2 and V3, tub being the B 
meson mass. £2 (£3) are the polarization vectors of the meson V% (V3). The amplitudes MS"^ corresponding 
to the polarization configurations a with both V2 and V3 being longitudinally polarized, and being transversely 
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polarized in the parallel and perpendicular directions are written as 

M° = {m%M L , m 2 B M N e* 2 (T)-et(T), -iM T e a ^ p e* 2a (T)e^(T)P 2j P 3p ) , 



(3) 



respectively. In the above expressions e(T) denote the transverse polarization vectors, and we have adopted the 
convention eoi23 = 1. 

Define the velocity t>2 = P-xlvny^ (v^ — Ps/mv 3 ) in terms of the V2 (V3) meson mass my 2 (m^). The helicity 
amplitudes, 



A L = -Gm 2 B M L , 
A» = G\f2m B M N , 



A± = Gm V2 mv 3 y/2[(v 2 ■ v 3 ) 2 - T\Mi 



with the normalization factor G = J G 2 F P C / (64irm B r), satisfy the relation, 

|A L | 2 + |^||| 2 + |A ± | 2 = 1. 
We also need to employ another equivalent set of helicity amplitudes, 

Mt 



H q = m B M L 



H± = m B M N T 



with the helicity summation, 



Y^M {a)t M (a) = \H \ 2 + \H+\ 2 + \H-\ 



The definitions in Eq. Q are related to those in Eq. © via 

G 



A, 



-GHq , 



-4 11 



V2 



(H+ + H-) 



A 1 



G_ 

V2 



(B+-H-) 



(4) 



(5) 



(6) 



(7) 



(8) 



The explicit expressions of the distribution amplitudes 



'pi 



, and (j) s p for a longitudinally polarized p meson, 



and 



b v p , and (j) a for a transversely polarized p meson are referred to |20t l2ll| . However, for the twist-3 



distribution amplitudes <fi p , (j) s p , <p v p , and cj) p , we adopt their asymptotic models as shown below: 

<p p (x) 

fT 

=3(2x- l) 2 



4> v P {x) 



^k=x{l-x) [l + 0.18C 2 3/2 (2:r-l) 



2y/2N c 



2^/2N t 

fT 



3f 

p x{l-x) l + 0.2C* 2 3/2 (2a;- 1) 



f2Nr. 



-[l + (2*-l) 2 ] 



2^/27V;4 

= (1 - 2x) 



3/ P 



Aj2N r 



(9) 
(10) 

(11) 
(12) 
(13) 
(14) 



with the decay constants f p = 200 MeV and /J = 160 MeV, and the Gegenbauer polynomial C^/ 2 {t) = 
3(5i 2 — l)/2. On one hand, the sum-rule derivation of light-cone meson distribution amplitudes suffer sizable 
theoretical uncertainty, so that the asymptotic models are acceptable. On the other hand, the asymptotic 
models for twist-3 distribution amplitudes were also adopted in QCDF 9], and the comparison of our results 
with theirs will be more consistent. 
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For the b — > d transition, the helicity amplitudes have the general expression, 

H h = v; b v ud + v; b v cd h<£> + v* b v td , (15) 

with h = or ±, and V's being the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. The amplitudes 
, , and are decomposed at LO into 

H ( u ] = m% (f p F e h +M h e + f B F> 1 + M h a ) , 

ntf = , 

= -m% (f P Fp h + M p e > h + f B F a p ' h + M p a ' h ) . (16) 
The LO PQCD factorization formulas for the B — ► pp helicity amplitudes associated with the final states p + p~ , 



P P H h 



F e h 4 (Ol) 

Me M h e4 (oi) 

^ -JT^ (02) 
■Ma (02) 



P + P~ 



(a't ) )+M h e Ja't ) ) 
nrF^(a^)+F^\a^) 



F P,h ph 

1 p. 1 e4 



F 



p,h 



M p a ' h M h a4 (af > + of ' + af > + + af ') + ^A^e 

°\ V2H ( h u) " 



(d) 

6 



P P 



Fi\ (01 + a 2 ) 
.M^K+c^) 
Fi 1 

Ma 



p+p" I ; 

' F e \ (a™ - a< d > + a<"> - af + a<"> - af) 

(af - 4 <d) + af ) - af > - af > + a^) + M% (a**> - o^) 



■Mr- ft i 

771 v^^"' 



^ F e h 4 (-02) 

A4* (-4) 

Jjrift (a 2 ) 

■M£ I -M^ 4 (4) 



pO I ..-;;"> 



^ 4 (4 d) )+^6(4 d) ) 

■Mf' h I -M^4 + a'} d) + + + 4 (<0 ) + (4 



(d) 

6 



(d) 
6 



TABLE III: LO B — > pp decay amplitudes with rfr = (1) for the longitudinal (transverse) components. 
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p + p°, and p°p° are summarized in Table Hill The Wilson coefficients for the factorizable contributions, and 
a '(i) f or th e nonfactorizable contributions can be found in Q , where q = u or d denotes the quark pair produced 
in the electroweak penguin. 

The explicit expressions of the LO factorizable amplitudes F® 4 a4 a6 and of the LO nonfactorizable amplitudes 
Me 4 e6 a4 a6 are similar to those for the B — > PP decays Q but with the replacements of the distribution 
amplitudes and the masses, 

<j> A (x)^cj>(x), 4 P (x)^ ^ (x) , f (x) - </>' (x) , 

mo2 -* -m p , m 03 -> m p . (17) 

In the above replacement mo2 (^03) is the chiral enhancement scale associated with the pseudo-scalar meson 
involved in the B —> P transition (emitted from the weak vertex), and m p = 0.77 GeV the p meson mass. 
Note that the amplitude F® e from the operators O5-S vanishes at LO. The LO factorization formulas for the 
transverse components are collected in Appendix A, whose relations to F^ and to in Table ITTT1 follow 
Eq. JBJ). For example, the amplitude F^ is given by 

F e4 = F Ne4 T ■ (18) 



B. NLO Corrections 



The vertex corrections to the B — ► pp decays modify the Wilson coefficients for the emission amplitudes in 
the standard definitions into 



ai(fx) 
a 2 (p) 
a.i(p) 



+ — ; — C F ——V 1 {p) , 



4tt 



02(W H -. C F — - — V 2 {p) , 



4tt 



Nr 



a s {p) Cj±x{p) 
a i(N + — i — c f — 77 Vi{p) , 2 = 3-10 



4tt N c 

where Vi(p) in the NDR scheme are in agreement with those in |22j | for the longitudinal component, 

_ is /"' 

p, 



Vt(p) 



121n^-18+^ 

fp 



dxcj)p(x)g(x) 



for i = 1 - 4, 9, 10 , 



, m b 2 y f 2N r 
-12 In— + 6- 



2y/2W c 

J a 



dx 4> p {x) g(l — x) , for i = 5, 7 



da;^"(a;) [-6 + fc(a:)] , 



for i = 6, < 



ip Jo 

and with those in |23j for the transverse components, 



, to,, 2^/2lV7 
12 In— -18 



for i = 1 - 4, 9, 10 , 



, m b „ 2y/2N r 

-12 In— + 6- 



dx[<f> v p (x)±cf> a Jx)}g(l-x) , forz = 5,7. 



P /p 
iate 

corrections introduce the additional contributions resulting from the penguin operators O^s, 



(19) 



(20) 



(21) 



We do not show V^ s , because of the associated factorizable emission amplitudes F^ = 0. Moreover, the vertex 





P 


f P F e PM " 


- fpF e P ' h 


p' 


-ffi: 


f p F P ' h - 


- f P F?< h 


p°p° : 


f P F e P ' h - 





) _ (d) 
6VC a 6VC I ' 



(22) 
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where the arguments etgvc represent only the vertex-correction piece in Eq. I|19l) . 

Taking into account the NLO contributions from the quark loops and from the magnetic penguin, the hclicity 
amplitudes are modified into 



(a) 



p p . H h -> # h + m^-M^ , H h - m B M h - to b M^ 

pV : i^"' c ' 4) -> iJ< u ' c,t) , (23) 



where -Ml , M.u , M$ , and denote the up-loop, charm-loop, QCD-penguin-loop, and magnetic-penguin 
corrections, respectively. The magnetic-penguin contribution to the B — + PV modes was computed in |2 1|. 

M^' Ctt) and M[ 9) for h = are similar to those for the B — ► PP decays Q with the replacements in Eq. 1|17|) . 
Those for the transverse components are presented in Appendix A. 

The choices of the B meson wave function, of the B meson lifetimes, and of the CKM matrix elements, 
including the allowed ranges of their variations, are the same as in Q . We vary the Gegenbauer coefficients in 
4> p and in <^>T by 100% as analyzing the theoretical uncertainty. The resultant B — > p form factors at maximal 
recoil, 

A = 0.311°$ , A 1 = 0.21±S;8| , V = 0.26±g;gg , (24) 

associated with the longitudinal, parallel, an d perpe ndicular components of the B ^ pp decays, respectively, are 
similar to those derived from QCD sum rules |2g,[2g, and almost the same as adopted in the QCDF analysis [27j |. 
Compared to [2f| , one-loop radiative corrections to the two-parton twist-3 contributions have been considered 
in |2(j . The central value of the form factor V in Eq. 12411 is a bit smaller than those in [2^, [2(| . We emphasize 
that this difference is not essential, since the perpendicular component corresponding to V contributes roughly 
less than 10% of the total B — ► pp branching ratios as shown below. 

The PQCD results for the B — > pp branching ratios, together with the BABAR and Belle data, are listed 
in Table HVI It is obvious that the NLO PQCD values are consistent with the data of the B a — > p T p^ and 
B^ — > p^p branching ratios. The color-suppressed tree amplitude is also enhanced by the vertex corrections 
here, but the ratio \C/T\ w 0.2 for the longitudinal component, similar to that in the B — > irir decays, is 
still small. However, the central value of the predicted B° — > p°p a branching ratio has almost saturated the 
experimental upper bound. We conclude that it is unlikely to accommodate the measured B° — ► 7r°7r°, p°p° 
branching ratios simultaneously in PQCD. 

Mode BABAR fl] Belle [1] LP LOnlowc +VC +QL +MP +NLO 

B°^>p f p k 30±4± 5 22.8 ±3.8^2.6 27.8 26.1 25.2 26.6 25.9 25.3+^ ( ( l"gj 

B ± -> p ± p° 17.2 ±2.5 ±2.8 31.7 ± 7-l±|;f 13.7 16.2 16.0 16.2 16.2 16.0+ 1 g-' 1 > ( ( + 

B° -> p°p° < 1.1 — 0.33 0.56 1.02 0.62 0.45 0.92 + _l ) 1 5 ° 6 [t° 



TABLE IV: B -» pp branching ratios from PQCD in the NDR scheme in units of 10~ 6 . 

We obtain the direct CP asymmetries A C p(B Q -► p^ p^) = -0.02 (-0.07), A CP (B ± -> p ± p°) = 0.00 (0.00), 
and A C p{B° -> p°p°) = 0.56 (0.80), where the values (in the parentheses) are from LO (NLO) PQCD. We have 
also computed the polarization fractions. The NLO corrections have a minor impact on the B° — > p^p^ and 
B — > p^p° decays: their longitudinal polarization contributions remain dominant, reaching 93% and 97%, 
respectively. However, the B° — > p°p° polarization fractions are sensitive to the NLO corrections as indicated 
in Table [V] where the average longitudinal, parallel, and perpendicular polarization fractions, /||, and 
respectively, are defined by 



B(B° -> p°p°) L , lU + B(B° -> pVh,i ± 
B(B° -» p°p°) + B(B° -» p°p°) 



s 



The av erag e longitudinal polarization fraction of the B° — ► p°p° decays was also found to be smaller in LO 
PQCD [181 llSj . It is easy to understand the changes due to the NLO effects. As stated before, the color- 
suppressed tree amplitude, being the main tree contribution in the B° — ► p°p° decay, is enhanced b y th e vertex 
corrections. The B° — ► p°p° polarization fractions should then approach the naive counting rules [271 I28L l2fll|: 
ft, ~ 1 and /|| ~ f± ~ A 2 obeyed by a tree-dominated decay, where A w 0.22 is the Wolfenstein parameter. 



Mode ft fj f ± 

B° -> p°p a 0.71 (0.67) 0.14 (0.15) 0.15 (0.18) 

B° -» p°p° 0.09 (0.79) 0.45 (0.10) 0.46 (0.11) 

Average 0.23 (0.78) 0.38 (0.11) 0.39 (0.11) 



TABLE V: LO and NLO (in the parentheses) polarization fractions of the B a — > p°p° decays from PQCD. 



III. JET FUNCTION IN SCET 

In this section we investigate the resolution to the B — > tttt puzzle claimed in QCDF with the input of the NLO 
jet function from SCET (lQ|. The leading-power SCET formalism has been derived for two-body nonleptonic B 
meson decays ^lj. However, there exist different opinions on the calculability of the hard coefficients in SCETn, 
one of which is the jet function characterized by a scale of O^m^K). In the jet function is regarded as being 
incalculable, and treated as a free parameter. Together with other hadronic parameters, it is determined by 
fitting the SCET formalism to the B — ► tttt data. Therefore, the large ratio \C/T\ obtained in ^1 is an indication 
of the data, instead of coming from an explicit evaluation of the amplitudes. In this analysis the QCD penguin 
amplitude, receiving a significant contribution from the long-distance charming penguin 30], was also found to 
be important. Similarly, the large charming penguin, as one of the fitting parameters in SCET, also arises from 
the data fitting. A global analysis of the B — > tttt, ttK decays based on the leading-power SCET parametrization 
has been performed recently in j^l), where a smaller branching ratio B(B° — > 7r 7r°) « 0.84 x 10~ 6 was obtained. 

A plausible mechanism in SCET for enhancing the ratio \C/T\ was provided in [ljj: f ne jet function could 
increase the nonfactorizable spectator contribution to the color-suppressed tree amplitude C at NLO. This 
significant effect was implemented into QCDF [10| . Because of the end-point singularities present in twist-3 
spectator amplitudes and in annihilation amplitudes, these contributions have to be parameterized in QCDF 
Different scenarios for choosing the free parameters, labelled by "default", "SI", "S2", ■ ■ ■. "S4", were proposed 
in [221 . As shown in Table IVT1 the large measured B° — > 7r°7r° branching ratio can be accommodated, when 
the parameter scenario S4 is adopted. It has been emphasized in the Introduction that the same formalism 
should be applied to other decay modes for a check, among which we focus on the quantities sensitive to C: the 
B — ► ttK direct CP asymmetries and the B° — ► p°p° branching ratio. 

The QCDF formulas for the B — * VV decays with the NLO contributions from the vertex corrections, the 
quark loops, and the magnetic penguin can be found in jUI^,!^, which appear as the 0(a s ) terms of the 
Wilson coefficients a,, i = 1, • • • , 10. The vertex corrections are the same as in Eqs. (|20|) and (|21ll . Note that 
the expressions of the Wilson coefficients a^g differ between [2?} and [23l l32l|: clq$ for both the longitudinal 
and transverse components in [2j|,|32 do not receive any 0(a s ) correction. We disagree on this result as shown 
in Eqs. (|20|l and i|23[l . Hence, we adopt the expressions in [27| for the contributions from the quark loops, 
the magnetic penguin, and the annihilation. We also employ the B — > p form factor values in 271. Since the 
spectator amplitudes were not shown explicitly in j2?|], we use those from |2^|. The parameter sets default and 
S4 have been defined for the B — ► PP decays [2^, but have not for the B — ► VV ones. Therefore, we assume 
that the parameters for the latter are the same as for the former in the following analysis. Fortunately, the 
predicted B — > p°p° branching ratio is insensitive to the variation of the annihilation phase (f>A, which is one of 
the most essential parameters in QCDF: varying <j>A between and 2-7T, the B° — > p°p° branching ratio changes 
by less than 10%. 

The jet function ju derived in is relevant to the B — > PP decays and to the B — » VV decays with 
longitudinally polarized final states. The jet function j± is relevant to the B — * VV decays with transversely 
polarized final states. These jet functions apply not only to the color-suppressed tree amplitudes, but to the 
color-allowed tree and penguin amplitudes, which are free of the end-point singularities. We mention that the 
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NLO corrections to the hard coefficients of SCET/ have been derived in [33L l34| . This new piece modifies the 
QCDF outcomes slightly, comparing the color-allowed and color-suppressed tree contributions obtained in [To| 
and in |33j | . Hence, we consider the NLO correction only from the jet function for simplicity. Furthermore, 
since the jet function enhances the color-suppressed tree amplitude, the B° — ► p°p° polarization fractions are 
expected to approach the naive counting rules. That is, the longitudinal component dominates. This tendency 
has been confirmed in PQCD as indicated by Table El To serve our purpose, it is enough to evaluate only the 
B — > plPl branching ratios here. 





Mode 


Data [1] 


default, LO jet default, NLO jet 


S4, LO jet 


S4, NLO jet 


B ± 




5.5 ±0.6 


6.02 (6.03) 


6.24 (6.28) 


5.07 (5.07) 


5.77 (5.87) 


B° 


— > ir^n^ 


4.9 ±0.4 


8.90 (8.86) 


8.69 (8.62) 


5.22 (5.17) 


4.68 (4.58) 


B° 


-» 7T°7r° 


1.45 ±0.29 


0.36 (0.35) 


0.40 (0.40) 


0.72 (0.70) 


1.07 (1.13) 


B ± 




24.1 ± 1.3 


20.50 (19.3) 


20.13 


21.60 (20.3) 


20.50 


B ± 




12.1 ±0.8 


11.79 (11.1) 


11.64 


12.48 (11.7) 


12.02 


B° 




18.9 ±0.7 


17.33 (16.3) 


17.21 


19.60 (18.4) 


19.23 


B° 




11.5 ± 1.0 


7.49 ( 7.0) 


7.41 


8.56 ( 8.0) 


8.36 


B ± 




19.1 ±3.5 


18.51 


19.48 


16.61 


18.64 


B° 


^pIpI 


25.2±|:f 


25.36 


24.42 


18.48 


16.76 


B° 


- pIpI 


< 1.1 


0.43 


0.66 


0.92 


1.73 



TABLE VI: Branching ratios from QCDF with the input of the SCET jet function in units of 10 . The values in the 
parentheses are quoted from [TH. l22j for comparison. The data for the B — > pp decays include all polarizations. 



Mode Data [1] default, LO jet default, NLO jet S4, LO jet S4, NLO jet 



B ± 




1±6 


-0.02 (-0.02) 


-0.02 


-0.02 


(- 


-0.02) 


-0.02 


B° 


— > TT^IT^ 


37 ± 10 


-6.57 (-6.5) 


-6.65 


10.60 


( 


10.3) 


10.91 


B° 


-» jrV 


90+40 


44.67 ( 45.1) 


41.95 


-19.58 


(- 


-19.0) 


-18.48 


B ± 


7r±A'° 


-2 ±4 


0.84 ( 0.9) 


0.85 


0.29 


( 


0.3) 


0.29 


B ± 




4±4 


6.88 ( 7.1) 


7.04 


-3.53 


(- 


-3.6) 


-4.08 


B° 




-10.8 ± 1.7 


4.28 ( 4.5) 


4.24 


-4.06 


(- 


-4.1) 


-3.89 


B° 




2 ±13 


-3.15 (-3.3) 


-3.37 


0.78 


( 


0.8) 


1.60 



TABLE VII: Direct CP asymmetries from QCDF with the input of the SCET jet function in percentage. The values in 
the parentheses are quoted from (22^1 for comparison. 



The predictions for the B — > tttt, ttK, plPl decays from QCDF with the input of the SCET jet function are 
summarized in Tables fVTI and IVTTl The values in the parentheses are quoted from ^(j f° r t ne B — ► 7T7t branching 
ratios, and from |22| for the B — ► hit direct CP asymmetries and for the B — ► ttK decays. The small differences 
between our results and those from [i"oll22| are attributed to the different choices of the CKM matrix elements, 
meson masses, etc. All the calculations performed in this work, except those of the B — > 7T7t branching ratios, are 
new. It is found that the scenario S4 plus the NLO jet function lead to the ratio C/T w 0.8, and accommodate 
at least the BABAR data of the B° — > 7r 7r° branching ratio. Nevertheless, the same configuration overshoots 
the experimental upper bound of the B° —> p°p° branching ratio apparently, implying that the color-suppressed 
tree amplitude is enhanced overmuch by the NLO jet function. Adopting the default scenario, QCDF satisfies 
the B° — > p°p° bound, but the predicted B° — > tt°tt q branching ratio becomes too small. We have surveyed 
the other scenarios, and found the results from SI and S3 (S2) similar to those from the default (S4). That 
is, it is also unlikely to accommodate the B — > tttt, pp data simultaneously in QCDF. The B — > ttK branching 
ratios are not affected by the NLO jet function, because the color-suppressed tree amplitude is still subleading 
in the penguin-dominated modes. The B^ — > p^p° L and B° — > branching ratios are not either, since they 

involve the larger color-allowed tree amplitude. 
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(a) (b) (c) 

FIG. 1: Contributions from inelastic FSI. 



Another indication against the resolution in 10] is given by the direct CP asymmetries of the B — > nK 
decays shown in Table IVTT1 both A C p(B ± -> tt ^) and A C p(B° -> ir^K^ deviate more from the data, a 
consequence expected from the discussion in As explained in Sec. II, the color-suppressed tree amplitude C 
needs to be roughly orthogonal to T" in order to have a vanishing Aqp{B ± — > 7r°i4T ± ). The NLO jet function, 
though increasing C, does not introduce a large strong phase relative to T". That is, C'/T' is large, but remains 
almost real as in the B — ► tttt case [T(j . This is exactly the same reason the B — > ttK puzzle can not be resolved 
in SCET [ll|3l|]: the leading-power SCET formalism demands a real ratio C'/T', such that a large C'/T' 
just pushes the SCET prediction for Acp(S ± — > 7r°if ± ), about —18% further away from the data. The 
direct CP asymmetry of the B° — > 71"°^° decays, whose tree contribution comes only from C", is sensitive to 
the NLO jet function as indicated in Table IVm The direct CP asymmetries of the B — > tttt decays, which are 
tree-dominated, are relatively insensitive to the NLO jet function. 



IV. DISCUSSION 



Before concluding this work, we comment on and compare the various analyses of the FSI effects in the 
B — > irK, tttt decays. The tiny B° — > 7r°7r° branching ratio obtained in perturbative calculations naturally 
leads to the conjecture that FSI may play an essential role. Though the estimate of FSI effects is very model- 
dependent, the simultaneous applications to different decay modes can still impose a constraint. The FSI effects 
from both the elastic and inelastic channels have been computed in the Regge model for the B — > tttt decays [l5| 
and for the B — > VV decays j^- The conclusion is that FSI improves the agreement between the theoretical 
predictions and the experimental data, but does not suffice to resolve the B — ► tttt puzzle: the B° — > 7r°7r° 
branching ratio is increased by FSI only up to 0.1-0.65 .15]. Moreover, the inelastic FSI through the long- 
distance charming penguin was found to be negligible in the B — > irn decays, though it might be important in 
the B — > irK ones. The reason is that the contribution from the intermediate DD states is CKM suppressed 
in the former compared to the D S D states in the latter. This observation differs from that in |lllll2| . where a 
significant charming-penguin contribution was claimed. We have pointed out in Sec. Ill that the large charming 
penguin in 0, 0] is a consequence of fitting the SCET parametrization to the data. 

The inelastic FSI has been also evaluated as the absorptive part of charmed meson loops shown in Figs. Wi'' 1 ) 
and^b) 14] . The two unknown cutoff parameters, appearing in the form factors associated with the three- 
meson vertices, were fixed by the measured B — ► ttK branching ratios. Note that these parameters should be 
the same for B — ► irK and B — > tttt in the SU(3) limit. Applying the same formalism to the latter, FSI can 
not resolve the B — »• tttt puzzle, even allowing reasonable SU(3) breaking effects for the cutoff parameters. This 
result is understandable: the absorptive amplitudes from Figs, da) and ^b) are more or less orthogonal to 
the short-distance QCD penguin amplitudes in the B — ► tttt decays, so that their effect is minor. Hence, the 
conclusion in 14] is the same as in |l5j |. That is, the charming penguin is not enough to explain the observed 
B — > tttt branching ratios. 

Then additional dispersive amplitudes must be taken into account in [l4j]. Those from Figs.^a) and^b), 
though calculable in the framework of [l4|. were not considered. If considered, they, also contributing to 
the B — > ttK decays, would change the earlier predictions. Therefore, a brand new mechanism, the dispersive 
amplitude from the meson annihilation DD — > tttt shown in Fig.^c), was introduced. There is no corresponding 
diagram for the B — > ttK decays. However, this amplitude is beyond the theoretical framework, i.e., it can not 
be expressed in terms of the Feynman rules derived in |14| . The four free parameters, namely, the two cutoff 
parameters involved in Figs.^a) and^b), and the real and imaginary contributions from Fig. ^c), were then 
determined by the four pieces of the B — > 7T7t data: the three branching ratios and the direct CP asymmetry 
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Acp(B° — > n^^). That is, the B° — ► 7r°7r° branching ratio has been treated as an input. The point of 0] is 
to predict the direct CP asymmetries of the B° — ► 7r°7r° and S ± — > vr^ decays, using the parameters fixed 
above. 

The rescattering among the final states of the B — > PP decays with P = n, K, and r\ has been studied in [T^j ], 
These elastic FSI effects were parameterized in terms of two strong phases, which, together with the B — > 7r 
and B ^ K form factors, and the chiral enhancement scale, were then determined by a global fit to the data, 
including the measured B° — > 7r°7r° branching ratio. Nevertheless, the feature of the elastic FSI effects, i.e., the 
correlated decrease and increase of the B° — ► tt^tt^ and B° — > 7r°7r° branching ratios, respectively, was noticed 
H3l. A FSI phase difference between the two B — > 7T7t isospin amplitudes with J = 0, 1 has been introduced 
in |33, which was then varied to fit the B — > 7T7t data. Therefore, no explanation for the large -B — > 7r 7r° 
branching ratio was provided from the viewpoint of FSI. 

There exist other global fits based on different parametrizations for the charmless B meson decays. For 
example, the large rati o C/ T was extracted by fitting the quark-amplitude parametrization to the B — > mr 
data [sH Issl Isfll l4Cl liH 1421 liM . No responsible mechanism was addressed, though the largeness of C was 
translated into the largeness of the QCD penguin with an internal t quark and/or of the exchange amplitudes 
in |39j ]. The QCDF formalism, in which the twist-3 spectator and annihilation amplitudes with the end-point 
sing ularities were parameterized as mentioned in Sec. Ill, has been implemented into a global fit to the data 
|45l 1461 |47|. To reach a better fit, the free parameters involved in QCDF must take different values for the 
B — > PP, PV, VP modes. These parameters have been tuned to account for the B — > irir data in 36]. As 
emphasized before, the analysis must be also applied to other modes in order to obtain a consistent picture: the 
parameters preferred in (3(| lead to a large real C/T, which is not favored by the data of the B — > nK direct 
CP asymmetries as stated in Sec. III. 

After carefully investigating the proposals available in the literature, we have found that none of them can 
really resolve the B — ► irir puzzle. The NLO PQCD analysis has confirmed that it is unlikely to accommodate 
the B — > 7T7T, pp data simultaneously (the NLO PQCD predictions are consistent with the B — > pp data). 
The B — > 7T7T decays have been studied in the framework of light-cone sum rules (LCSR) 01 1 where a small 
B° — > 7r°7r° branching ratio was also observed. Since there is only little difference between the sum rules for 
the B — > 7T7T and B — > pp modes, we expect that the conclusion from LCSR will be the same as from PQCD. 
The resolution with the input of the NLO SCET jet function into QCDF [l0( does not survive the constraint 
from the B — ► pp data, and renders the B° — ► n^K^ and B ± — > tt^K^ direct CP asymmetries deviate more 
away from the measured values. We conclude that the B —> pp data have seriously constrained the possibility 
of resolving the B — * tttt puzzle in the available theoretical approaches. 
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supported by the National Science Council of R.O.C. under Grant No. NSC-94-2112-M-001-001, by the Taipei 
branch of the National Center for Theoretical Sciences, and by the U.S. Department of Energy under Grant 
No. DE-FG02-90ER40542. 

APPENDIX A: TRANSVERSE HELICITY AMPLITUDES 

In this Appendix we present the factorization formulas for the transverse helicity amplitudes: 

Fjve4(a) = l&irCFTn 2 B I dx\dx% I b\db\ bidb% 4>b{%i, &i) 
Jq Jq 

xr 3 I [^(ZJ) + 2r2$(Z2) + r 2 x 2 (^(Zj) + $(25))] E e {t)h e {A,B,b 1 ,b 2 ,x 2 ) 

+r 2 (^ 2 (x^)-^(x^))E e (t')h e (A',B',b 2 ,b 1 ,x 1 )^ , (Al) 

/>1 />oo 

Frei{o) = 32irCFm% / dx\dx 2 I b\db\ b 2 db 2 4>b{x\, b\) 
Jq Jq 
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x r 3 | [4>l (55) - 2r 2 05 (55) - r 2 x 2 (05 (55) + 05 (35) )] E e (t ) h e ( A, B , h , b 2 , x 2 ) 



+r 2 ( ( p v 2 (xZ)- ( j>z(xZ))E e (t')h e (A',B',b2,b 1 ,xi)j , (A2) 
F Ne6 {a) = F Te6 (a) = 0, (A3) 

FNai(a) = l6TrCFm 2 B r 2 r 3 / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

X I [(1 - X 3 ) ($(3a)$(Z3) + 05(35)05(35)) + (1 + *s) (05(35)05(35) + 05(35)05(35))] 
xE a (t)h e (A,B, b 2 ,b 3 ,x 3 ) 

- [(2 - x 2 ) (05(35)05(35) + $(282)05(33)) - *2 (05(35)05(23) + 05(35)05(35))] 



xE a (t')h e (A',B',b 3 ,b 2 ,x 2 )^ , 

WirCpin^ r 2 r 3 / dx 2 dx 3 / 62^2 &3<#>3 
Jo Jo 



x I [a* (05(5^) - 05(55)) (05(35) - ^(^)) + (05(35) + 05(35)) (05(3*) + 05(35))] 
xE a (t)h e (A,B,b 2 ,b 3 ,x 3 ) 

- [(1 - x 2 ) (05(55) + 05(55)) (05(55) + 05(35)) + (05(35) - 05(35)) (05(35) - 05(35))] 
xE a (t')h e (A',B',b 3 ,b 2 ,x 2 )} , (A4) 

r-l />00 

FTai{a) = - 32irC F m 2 B r 2 r 3 / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x I [(1 - x 3 ) (05(55)05(35) + 05(35)05(35)) + (1 + x 3 ) (05(55)05(35) + 05(35)05(35))] 
xE a (t)h e (A,B,b 2 ,b 3 ,x 3 ) 

- [(2 - a; 2 ) (05(55)05(35) + 05(35)05(35)) - x 2 (05(55)05(35) + 05(35)05(35))] 
xE a {t')h e (A',B',b 3 ,b 2 ,x 2 )^ , 

/>1 poo 

= Z2irCFm B r 2 r 3 I dx 2 dx 3 I b 2 db 2 b 3 db 3 
Jo Jo 

x I [x 3 (05(55) - 05(35)) (05(35) - 05(35)) - (05(35) + 05(35)) (05(35) + 05(35))] 
xE a (t)h e (A, B,b 2 , b 3 ,x 3 ) 

+ [(1 - x 2 ) (05(55) + 05(35)) (05(35) + 05(35)) - (05(35) - 05(35)) (05(35) - 05(35))] 
xE a {t')h e (A',B',b 3 ,b 2 ,x 2 )} , (A5) 

F^ a %{o) = 32irCFm 2 B / dx 2 dx 3 I b 2 db 2 b 3 db 3 
Jo Jo 

x |r 2 (05(55) + 05(35)) ^(x^)E a (t)h e (A, B, b 2 ,b 3 , x 3 ) 

+r 3 ^(xZ) (05(55)- 05(55)) E a (t')h e (A',B', b 3 ,b 2 ,x 2 )} , (A6) 
F Ta6 (a) = 2F Na6 (a) , (A7) 



^/2N~ C f 1 f°° 

2,2-kCf — -rr-^-m 2 B r 3 I dx\dx 2 dx 3 I b\db\b 3 db 3 (j)B{x\,bi) 

Jo Jo 

x{(i-x 3 )<fi(iB2) (M^)~^))K(t)K(A,BMM) 

+ [x 3 <j>l{xl) ($(3s) - 0§(2s)) - 2r 2 (x 2 + x 3 ) (^(^)^(^) + ^(^)^(^))] 
xE' e {t')h n {A',B'MM)} , 

32itCf — —^m B r 3 I dx\dx 2 dx 3 \ b\db\b 3 db 3 §B{x\,b\) 

Jo Jo 

xj(i- 33)^(22) {4>im)~^))E' e (t)h n {A,BMM) 

+ [x^Hxi) (^(33) - 0§(33)) - r 2 (x 2 + a*) (^(3a) + $(3a)) (0§(S3) + $(3*)) 
-r 2 (x 2 + x 3 ) (0^) - <%(X2)) ($(3*) - 0§(53))] ^(t'K^'.S', 61,63)} , 

MttCf — TT-^-m%r3 / dx\dx 2 dx 3 I bidbi b 3 db 3 0b (zi, 61) 
Jo Jo 

xj(l-z 3 )0;r(3i) (0§(S3)-05(lB3))^ e (t)/l n (A,B ) 61,63) 

+ [.T 3 2 n (3i) (0§(33) - 05(33")) + 2r 2 (z 2 + x 3 ) (0^(3^)0^(3^) + 0^(3^)0^(3^))] 

x^(t')/i„(^,B',6i,6 3 )} , 

QA-kCf — T—^-m B r 3 I dx\dx 2 dx 3 I b\db\ b 3 db 3 4>b{x\, 61) 
Jo jo 

x|(l- ^3)0^ (3i) (0§(S3) - 05(33-)) ^(^(AB, 61, 63) 

+ [x 3 cg(x^) (^[(33) - 05(33")) + r 2 (x 2 + a*) (0^(3i) + $(3a)) (0§(S3) + 05(3^)) 

-r 2 (x 2 + 13) {4>l{xi) - 2 , (3 2 ")) ($(3-3) - 05(33"))] ££(i')M^'> B', 61, 63) 

-M N ei{a') , 
-M T ei{a!) , 

— 32-kCf — T-f—rri B r 2 I dx\dx 2 dx 3 \ b\db\ b 3 db 3 (/>b(xi, 61) 

N c Jo Jo 

xx 2 (0^(3^) + $(3*)) <fl{xt) [E' e (t)h n (A, B, 61, 63) + ^(i')M^ B', 61,63)] , 

2/WjVe6(a') , 

-64ttCf at C m 2 B r 2 r 3 I dxidx 2 dx 3 I bidbib 3 db 3 4> B (xi,bi) 

N c Jo JO 

x (0^(^)0^(23) + <P%(x 2 -)r 3 (M))E' a (t')h n (A',B', 63, 61) , 

— 32ttCf — TT-^-rn 2 B r 2 r 3 / dx\dx 2 dx 3 I 61 d6i 63(^63 0s(a;i, 61) 

Jo Jo 

x( (0^) + 02(32")) ($(3*) + 05(3i)) + (0^(3^) - 0^)) ($(3*) - 05(31)) 



x^OMA'.B', 6 3 ,&i) 
/27N 



1287tCf Ar c m 2 B r 2 r 3 / dx\dx 2 dx 3 I b\dbib 3 db 3 <f>B{x\,bi) 
Jo Jo 
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x (^(^)^(^) + ^(^)^(^))^(f')/in(A', J B / ,6 3 ,bi) , 

MttCf at c m B r 2 r 3 / dx\dx 2 dx 3 / &i<#>i b 3 db 3 <Pb(xi, h) 
Jo Jo 

x| (^(^) + 0S(32)) (^3(23) + $(23)) - (^(35) - ^(22)) (0S(33) - 03(^3)) 

x£:(i')^(^'B>3,ti) , 
A^Ar a4 (a') , 

Al T a4(a') , 



M T ab{a') 

MNa&{a) = 32itCf V „ C m 2 B I dx\dx 2 dx 3 I bidbib 3 db 3 4>B{xi,bi) 

Jo JO 

x j [r 3 x 3 <fi(x2) (05(3*) - $(3*)) - r 2 (l - x 2 ) (^(^) + ^(^)) ^(a*)] 
x££(t)M4-B,&3,&i) 

+ [r 3 (2 - x^fe) (05(3*) - $(3*)) - r 2 (l + x 2 ) ($(35) + ^(^)) ^(a*)] 
x£:(t')ft n (i',5',6 3! 6i)} , 
A^Tae(a') = 2Al7Va6(a') • 



(A15) 
(A16) 
(A17) 



(A18) 
(A19) 



The quark-loop corrections M^ T for q = u, c, and t, and the magnetic-penguin corrections M^ T to the 
transverse components are written as 

C 2 



M 



(?) 



C 2 f 1 f°° 
= — l6m 2 B F r 3 / dx\dx 2 dx 3 \ b\db\b 2 db 2 $ B {x\,b\) 
\/2N c J J 

x{ [^m {<t>im) +$(33)) 

+r 2 (2 + 3:2)^(^)^(53) + r2a:2^(^)0S(23) + ^2^(^2)^(53)] 
xE^(t q ,l 2 )h e (A,B,b u b 2 ,x 2 ) 
+ r 2 4>l{xi]<j> v 3 {x^E^{t' q , l' 2 )h e (A',B', b 2 ,h, Xl )}, 



M 



(9) 



, 



(A20) 
(A21) 



(9) 

AT 



16m 



/2iv: 



dx\dx 2 dx 3 



c JO 



M 



(s) 



/•OO 

/ bidbib 2 db 2 b 3 db 3 4> B {xi,bi) 
Jo 

{ [-r 2 (l - x 2 ) (^(33) + ^(33)) ^(283) - r 3 (l + x 2 )x 3 cg(x^) (05(33) - ^(33)) 
-r 2 r 3 (l - z 2 ) ($(3*) + $(3*)) (05(^3) + 05(3*)) 

-r 2 r 3 .T 3 (l - 2x 2 ) (^(3i) - ^(Za)) (05 (Z3) - <A 3 (23~))] E g (t q )h g (A, B, C, bi,b 2 , b 3 , x 2 ) 

~r 2 r 3 x 3 {<&{X2) - <(> a 2 (x2)) (05(3*) - $(3*)) E g (t' g )h g (A' , B' , C , b 2 , h, b 3 , a*)} , (A22) 

C 2 f 1 f°° 

F / dx\dx 2 dx 3 I b\db\ b 2 db 2 b 3 db 3 4>b(xi, b\) 
/2N C Jo Jo 

x { [-r 2 (l - x 2 ) (05(3*) + 3 l (3*)) 0^(3*) - r 3 (l + x 2 )x 3 <g{x^) (05(3*) - 3 t (3*)) 

+r 2 r 3 (l - x 2 ) (4%(X2) + 02(3*)) (0g(3*) + 3 s (3*)) 

-r 2 r 3 x 3 (l - 2x 2 ) {4> v 2 {xZ) - <f>%(xl)) (05(3*) - 3 (3*))] E g (t q )h g (A, B, C, b\,b 2 , b 3 , x 2 ) 
-r 2 [(1 - x 2 )(%(x2)<t%(x3) - r 3 (l - 2x2)^(3^)^(33") - r 3 $(3*)$(3j3) 

-r 3 ^ (32)05(3s) + r 3 05 (3^)0^(3^)] B fl (t , ,)/i fl (^', £?', C", 6a, 6 1} b 3 , a*)} . (A23) 
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